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C/3 ■ ABSTRACT 

Aims. We aim to study in detail the peculiar mineralogy and structure of the circumstellar environment of two binary post-AGB stars, 
EPLyr and HD 52961. Both stars were selected from a larger sample of evolved disc sources observed with Spitzer and show unique 
solid-state and gas features in their infrared spectra. Moreover, they show a very small infrared excess in comparison with the other 
sample stars. 

Methods. The different dust and gas species are identified on the basis of high-resolution Spitzer-IRS spectra. We fit the full spectrum 
to constrain grain sizes and temperature distributions in the discs. This, combined with our broad-band spectral energy distribution 
and interferometric measurements, allows us to study the physical structure of the disc, using a self-consistent 2D radiative-transfer 
disc model. 

Results. We find that both stars have strong emission features due to CO2 gas, dominated by 12 C 16 02, but with clear 13 C I6 02 and even 
16 12 C 18 isotopic signatures. Crystalline silicates are apparent in both sources but proved very hard to model. EPLyr also shows 
evidence of mixed chemistry, with emission features of the rare class-C PAHs. Whether these PAHs reside in the oxygen-rich disc or 
in a carbon-rich outflow is still unclear. With the strongly processed silicates, the mixed chemistry and the low 12 C/ 13 C ratio, EPLyr 
resembles some silicate J-type stars, although the depleted photosphere makes nucleosynthetic signatures difficult to probe. We find 
that the disc environment of both sources is, to a first approximation, well modelled with a passive disc, but additional physics such 
as grain settling, radial dust distributions, and an outflow component must be included to explain the details of the observed spectral 
< ^* ) energy distributions in both stars. 
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1. Introduction The condensates in C-rich outflows show features of carbon- 
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s_| ' . r , r . ^„ r . , , species such as SiC, MgS or polycyclic aromatic hydrocarbons 

03 ■ The infrared spectra of post-AGB stars are often characterised by (pAHfj) ( | Ho n V et afll2oM E002; Peeter7IHD|2002). They 

strong spectral signatures. These are formed in the gas and dust- , , 7 ■ 7T <~ t <, ni 

. , ° F „ fe . ,. , . 6 r , are also characterised by an often very strong feature at 2 1 urn 

rich circumstellar environment (CE), which is a remnant of the (Kwok gt aJ mg . Volk gt aJ 1 999I | H onv et al.lr20 03 ) . The pho- 

strong mass loss that occurred during the previous asymptotic t heres of these 21 sources sh o w strong enha ncements of 

giant branch (AGB) evolutionary phase. The chemistry in this , , , „ . , 17; • . , |hnn/il hnmk 

6 . • • , , -,i carbon and s-process elements (e.g. Revmers et al. 2004, 2007) 

circumstellar environment is found to be oxygen-rich or carbon- and ^ 2 j star s are recogn ised as post-AGB carbon stars 

rich, depending on whether oxygen or carbon is more abundant. ( Van winckd & Reyniers 2000) 
The less abundant of the two will be locked in the very stable 

CO molecule that forms in the stellar photosphere. Some eV olved objects show, however, features of both O- 

Typical post-AGB outflow sources that have O-nch CE not rich and c . rich dust species in their spectra ^ called 

only show the well-known 9.7 and 18 Aim features of amor- j^^ed chemistry sources. This chemistry is detected in sev- 

phous silicates but also narrower features, arising fromcrxs- eral sources in a wide range of differ ent evolutionary stages, 

talline silicates, (e.g. I Waters et al. | [1996t | Molster et alj |2002|). Some examp i es inc i ude Herbig Ae stars, or AG B stars, such as 

I _ j ; : , ... , . . , . J-type carbon stars with silicate dust emission (Little-Marenin 

Based on observations made with the 1.2m Flemish Mercator tele- rr^ m |i ,„»^ i , . — ~ 7^ 

scope at Roque de los Muchachos, Spain, the 1.2m Swiss Euler tele- 1986 ' Moyd Evans| | 1990 ^. Others are red giants, for example 

scope at La Silla, Chile and on observations made with the Spitzer HD 233517 an evolved O-nch red giant with PAHs in a circum- 

Space Telescope (program id 3274), which is operated by the Jet stellar disc (|Jura et al.||2006|). Other examples are planetary neb- 

Propulsion Laboratory, California Institute of Technology under a con- ulae (PNe) with evidence of silicates and PAHs dKemper et al.l 

tract with NASA. ** NASA Postdoctoral Fellow " 120021: iGutenkunst et alj|2008l) . or the hydrogen-poor [WC] cen- 
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tral stars of PNe dWaters etal.lll998t ICohen et al.lll999l) . Also 
some M supergiants are associated wit h emission due to PAHs 
dSvlvester et alJll998tlsioan et al.ll2008l) . 

Post-AGB stars with evidence of mixed chemistry in- 
clude HD 44179, the central star of the carbon-rich Red 
Rectangle nebula dCohen et al.lll975l) . The central star is a bi- 
nary surrounded by a Keplerian O-rich circumbinary disc (e.g. 
Van Winckel et al.ll 19951: 1 Waters et alJl998tlMen'shchikov et all 
2002; Bui arrabal et alj|2005b . Here the formation of the disc is 
belie ved to have antedated the C-rich tra nsition of the central star 
(e.g. lCohen etai1l2004HWitt et al.ll2008l) . 

Studies have shown that these evolved binaries with cir- 
cumbinary discs are much more ab u ndant than anticipated 
dDe Ruyter et all 120061: [VanWinckell 120071) . Interferometric 
studies dPeroo et alj2006 , 2007ab prove that the discs are indeed 
very compact, with radii around 50 AU in the N-band. The discs 
are also the natural environment of the obs erved photospheric 
chemical depletion patter n in these stars dVan Winckel et all 
1998; Giridhar et al.l2000l) . due to chemical f ractionation by dust 
formation in the circumstellar environment dWaters et al.1 11992) 
and subsequent accretion of the gas component. The presence of 
a long-lived stable reservoir of dust grains also could allow for 
the observed strong processing of the silicate dust grains, both in 
size a s well as in crystallinity (Molster et al. 2002t iGielen et al] 
l2008h . 

Dusty RVTauri stars are a distinct class in the post-AGB 
stars. They cross the instability strip, and are therefore pulsatin g 
stars dJuralll986t Uura & Kahanelll999l: iDeRuvter et all 120051) . 
RV Tauri stars show large-amplitude photometric variations with 
alternating deep and shallow minima. The members are located 
in the high-luminosity end of the population II instability strip, 
and the photometric variations are interpreted as being due to 
radial pulsatio ns. Circum stellar dust emission was observed in 
many of them dJural 19861) . and this was generally acknowledged 
to be a decisive character to place these stars in the post-AGB 
phase of evolution. The grains in almost all dusty RV Tauri stars 
are, however, not freely expanding but likely also trapped in a 
disc dVan Winckel et alj|l999t|Pe Ruyter et al1l2005M2006l) . 

In this paper we focus on two peculiar post-AGB stars 
with RV Tauri pulsational characteristics: EP Lyr and HD 5296 1 . 
These stars show unique spectral features and have very small 
infrared excesses in comparison to the larger sample. 

The outline of the paper is as follows: We start with a short 
description of the programme stars in Sect. [2] In Sect. [3] we 
give an overview of the different observations and reduction 
strategies. The analysis based on the Spitzer spectra is given in 
Sect. [4] and subdivided in different subsections. Sect. 14.21 con- 
tains a description of the silicate dust features and the modelling 
of the Spitzer-IRS spectra. The CO2 gas features are discussed in 
Sect. l4.3l and the observed PAH features in EP Lyr in Sect. 14.41 In 
Sect.[5]we model the observed SEDs using a passive disc model, 
also constrained with MIDI interferometric measurements. The 
discussion of our different results and our conclusions are pre- 
sented in Sect. [6] 



2. Programme stars 

In our previous study we described and modelled the Spitzer- 
IRS spectra of 21 sources and found that the dust around these 
stars is all O-rich and on average highly crystalline (Giele rTet al.l 
2008). The two stars discussed here have the lowest Lir/L», re- 
spectively 12% and 3%, in the larger Spitzer sample, where an 
average of about 50% was found. The large infrared luminosity 
can be explained with a passive disc model, provided that the 



inner rim is clo se to the star and th e scale height of the disc is 
significant (e.g. iDeroo et al.lr2007al) . The low observed Lir/L» 
values of both stars point to a small disc scale height and/or a 
much larger inner gap, as it is unlikely that a disc is optically thin 
in the radial direction. Not only do they have the lowest Lir/L» 
values, both stars show unique spectral signatures in comparison 
to the larger sample. We therefore selected these objects for a 
more detailed analysis. 

2.1. EPLyr 



ISchnellerl d 193 lb discovered the variability of EP Lyr and clas- 
sified it as an RVb star. RVb stars are objects with a vari- 
able mean magnitude, in th e General Catalog ue of Variable 
Stars dKholopov et al.1 1 19991) . Other studies IZsoldosI 11995b 
iGonzalez et al.lfl997h classify the light curve as an RVa photo- 
metric variable, having a constant mea n ma gnitude, with a pe- 
riod of P — 83.46 days. Preston et al ] d 1963b class ify it as an 
RVB spectroscopic variable. Gonzal ez et alj d 1997b performed 
an abundance analysis on EPLyr where they deduced stellar 
parameters (see Table Q} and found the star to be metal-poor, 
oxygen-rich and severely depleted. Using the molecular lines 
found in the spectra, they also quantified the 12 C/ 13 C ratio to be 
9 ± 1 . In the radial velocity data there is also evidence that EP Lyr 
must have a stellar mass companion, but additional observations 
are necessary to determine the orbit. 

2.2. HD 52961 

HP 52961 is an RV Taur i like object, similar to class RVb objects 
dWaelkens et al] [199 1b). with a p hotometric variability of 72 
days due to clear radial pulsations dWaelkens et a l. 1991b) . The 
binarit y of HP 52961 was first reported bv | Van Winckel et al.l 
d 19951) and furthe r refined in I Van Winckel et al] d 1999b and 
iPeroo et al] d2006h . where an orbital period of P orh = 1297 + 7 
days and an eccentricity of e — 0.22 + 0.05 was found. On top 
of the stable photometric variation due to the pulsation, another 
long-term photome tric variation was detect ed, correlated with 
the orbital period. IVan Winckel et al] d 1999b conclude that this 
can be understood as caused by variable circumstellar extinction 
during the orbital motion. 

The star is a highl y metal-poor object with [Fe/H] = -4.8 
dWaelkens et al.lll991ab and has an extremely h igh zinc to iron 
ratio of [Zn/Fe] = +3.1 dVan Winckel et al.ll992b . The star is one 
of the most extremely depleted objects known. 

HP 52961 has been studied with mid-IR long-baseline inter - 
ferometry using the VLTI/MIPI instrument dPeroo et al] |2006). 
They find that the dust emission originates from a very small but 
resolved region, estimated to be ~ 35 mas at 8 fim and ~ 55 mas 
at 13y«m, likely trapped in a stable disc. The dust distribution 
through the disc is not homogeneous: the crystallinity is higher 
in the hotter inner region. 



3. Observations 

High- and low -resolution spectra of 21 post-AG B stars were 
obtain ed using the Infrared Spectrograph (IRS; iHouck et al] 
2004) aboard the Spitzer Space Telescope (We rner et al] [2004) 
in February 2005. The spectra were observed using combi- 
nations of the short-low (SL), short-high (SH) and long-high 
(LH) modules. SL (/i=5.3 - 14.5 /mi) spectra have a resolving 
power of R=A/ A A ~ 100, SH (,1=10.0 - 19.5 //m) and LH 
(,1=19.3 - 37.0/im) spectra have a resolving power of ~ 600. 
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Name 


a (J2000) 


6 (J2000) 


^eff 


logg 


[Fe/H] 


Porb 


e E(B - V0,o, 




d 






(h m s) 


(° ' ") 


(K) 


(cgs) 




(days) 




(%) 


(kpc) 




EPLyr 


19 18 17.5 


+27 50 38 


7000 


2.0 


-1.5 




0.51 ±0.01 


3 + 


3.2 + 


1.0 


HD 52961 


07 03 39.6 


+ 10 46 13 


6000 


0.5 


-4.8 


1310 


0.21 0.06 + 0.02 


12 ± 1 


1.6 + 


0.5 



Note: Listed are the name, equatorial coordinates a and d(J2000), effectiv e temperature r c ff, surface gravity logg, and metallicity [Fe/H] of our 
sa mple stars. For the model parameters we refer to lDe Ruvter et alj J2006I) . Also given are the orbital period and the eccentricity (see references 
in lDe Ruvter et alJl2006l : lGielen et al1l2007h . The total reddening E(B - V) t 0l , the energy ratio Ltr/L, and the calculated distance, assuming a 
luminosity of L, = 3000 ± 2000 L Q . 



Exposure times were chosen to achieve an S/N ratio of around 
400 for the high-resolution modes, which we complemented 
with short exposures in low-resolution mode with an S/N ra- 
tio around 100, using the first generation of the exposure time 
calculator of the call for proposals. 

The spectra were extracted from the SSC data pipeline ver- 
sion S13.2.0 products, using the c2d Interactive Analysis reduc- 
tion s oftware package (Kess ler-Silacci et ai]|2006t ILahuis et al.l 
120061) . This data processing includes bad-pixel correction, ex- 
traction, defringing and order matching. To match the different 
orders, we applied small scaling corrections. 



4. Spectral analysis 

4.1. General 

Alookatthe spectra of EP Lyr and HD52961 (Fig.[TJi show very 
rich spectra with quite different continuum slopes. EPLyr shows 
strong emission features at longer wavelengths, with peak emis- 
sion in the 20 /mi region, whereas HD 52961 is characterised by 
a strong 10 /mi emission feature on top of a much steeper con- 
tinuum. 

Common dust species found in oxygen-rich post-AGB 
stars are amorphous silicates, namely olivine and pyroxene. 
Amorphous olivine (Mg2 A Fe2(i- A )Si04, where < x < 1 
denotes the magnesium content) has very prominent broad 
features around 9.8/mi and 18/mi. Amorphous pyroxene 
(Mg v Fei_ A Si03) shows a 10/mi feature similar to that of amor- 
phous olivine, but shifted towards shorter wavelengths. Also the 
shape of the 18 /mi feature is slightly different. For EPLyr it 
is unclear whether there is a significant contribution of amor- 
phous silicates. Small amorphous silicates could contribute to 
the observed strong emission bump at 20 //m in EPLyr, but 
as there does not seem to be a 10 //m amorphous feature, the 
20 //m bump could be purely continuum dominated. HD 52961 
has clear strong emission of amorphous silicates at 10 /mi, but 
the profile shows complex narrow subfeatures. Very little contri- 
bution at 20 /mi is seen. 

Both stars show strong narrow emission features which can 
be identified as being due to crystalline silicates. The Mg- 
rich end members of crystalline olivine and pyroxene, forsterite 
(Mg2SiC<4) and enstatite (MgSiC^), show strong but narrow fea- 
tures at distinct wavelengths around 1 1.3- 16.2- 19.7-23.7-28 
and 33.6//m. For EPLyr the silicate emission only clearly starts 
longward of 18/im, where strong emission features around 19 — 
23 - 27 and 33 //m can be seen. HD 52961 has strong narrow 
features at 9.8 - 1 1 .3 /mi and a remarkably strong 16 /mi feature. 
If this strong 16//m band is only due to forsterite it has shifted 
considerably to shorter wavelengths. A significant 16 /mi feature 



EPLyr 




15 20 25 

Wavelength [/im] 



HD52961 



40 




15 20 25 30 35 40 
Wavelength [/im] 

Fig. 1. The combined Spitzer-IRS high- and low-resolution spec- 
trum of EPLyr (top) and HD 52961 (bottom). In grey we over- 
plot the estimated errors. 



is seen in seve ral evolved disc so urces but it is never as strong as 
in HD 52961 dGielen et alj|2008l) . 

EPLyr shows evidence of the presence of carbon-rich dust 
species with probable PAH identifications at 8.1 and 11. 3 /mi. 
The detection of PAH emission together with silicates is surpris- 
ing and only observed in a few other post-AGB sources. The 
analysis of the PAH features is given in Sect. 14.41 

The spectrum of EP Lyr shows a strong resemblance to that 
of IRAS 09245-6040 (Fig. 0, a silicate J-type carbon AGB 
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star dMolster et al.ll2001l:lGarcfa-Hernandez et al.ll2006l) . Silicate 
J-type carbon stars have surprisingly low 12 C/ 13 C ratios and 
do not show the typical s-process o verabundances seen in N- 
type carbon stars ( Abia & Isern 2000). The infrared spectrum of 
these stars shows features of both carbon- and oxygen-rich dust 
species. 

Of the silicate J-type carbon stars, onl y 10% show emis- 
sion bands due to c rystalline material dLlovd Evansl I199U 
Ohnaka & Tsuii|[l999l) . The formation history of these stars is 
still unclear, but the most promising scenario for the presence 
of silicates in these stars, is that they a r e binaries with an unde - 
tected companion dLlovd Evansl 1 1990t lYam amura et al. 2000). 
A disc is supposed to be formed when the primary was still an 
O-rich giant. After that the star underwent thermal pulses and 
evolved into a carbon star. The silicate disc could be either cap- 
tured from the wind dMastrodemos & Morris! [19991) or the re- 
sult of a phase of strong binary interaction in a narrow system. 
Systems with strong crystalline features i n their spectra, suc h 
as IRAS 09425-6040 or IRAS 18006-3213 (iDeroo et al.ll2007bl) . 
would then be a result of mass-transfer into a circumbinary sys- 
tem, whereas sources dominated by amorphous sil icates, such as 
V778Cygni dYamamura et alj|2000l) or BH Gem dOhnaka et alj 
2008), consist of a wide binary with a circumcompanion disc. 
To date, no orbits are known, howeve r, and direct evidenc e of 
binarity is found in a few objects only (Izum iura et alf 2008). 

For IRAS 09245-6040, the 12 C/ 13 C ratio is calculated to be 
15+6 dGarcfa-Hernandez et al.1 120061). In the ISO-SWS spec- 
trum features of C2H-2, HCN, CO, C3 and SiC are seen short- 
ward of 15 jt/m; after 15 fim the spectrum is domin ated by strong 
emiss ion features of Mg-rich crystalline silicates dMolster et alJ 
12001b . As in EPLyr, there is no evidence of a strong contribution 
of amorphous silicates. 

Finally in both EPLyr and HD 52961, clear CO2 gas emis- 
sion features are detected the 13-18 jum region. This is discussed 
in Sect.[ 



EPLyr 



4.2. Silicate dust emission 



We op timised the fitting procedure as discussed in iGielen et al.l 
(2008) for these two outliers, where we modelled the full Spitzer 
sample, consisting of 21 stars. In short we assume the flux to be 
originating from an optically thin region, so we can make linear 
combinations of the absorption profiles to calculate the model 
spectrum. In our previous modelling we found that, on average 
for the full sample, the best fit was obtained using relatively 
large grains (> 2 /im) in an irregular Gaussian Random Fields 
(GRF) dust model. For EPLyr and HD 52961 however, we al- 
ready found that using smaller grain sizes (< 2 /mi) improved 
the fit considerably. 

So we repeated the analysis for EP Lyr and HD 5296 1 , allow- 
ing for different dust shapes, grain sizes and Mg/Fe content in the 
amorphous grains. We tested Mie, GRF and DHS (Distribution 
of Hollow Spheres) dust models in grain sizes ranging from 0. 1 
to 4.0 yum. In order to test for the presence of Fe-poor amorphous 
dust, we perform the modelling both with pure Mg-rich amor- 
phous silicates (x — 1) and with the more standard Mg-Fe amor- 
phous silicate dust (x - 0.5). For a detailed description of the 
fitting routine we refer to Gielen et al. (2008). The results of the 
fitting can be found in Table [2] As for EPLyr the silicate signa- 
tures only appear after 18 //m; we only fit this part of the Spitzer 
spectrum. 

The x 2 values of our fitting (Table [2]l are still qu i te hig h for 
HD 52961 but, confirming the result of lGielen etakl d2008l) . we 
can already tell that for both stars Mie theory is not a good dust 



0.5 



0.4 



0.3 



0.2 




25 30 
A (Aim) 

HD52961 




Fig. 2. Best fits for EPLyr and HD 52961. The observed spec- 
trum (black curve) is plotted together with the best model fit (red 
curve) and the continuum (black solid line). Forsterite is plotted 
in dash-dot lines (green) and enstatite in dash-dot-dotted lines 
(blue). Small amorphous grains are plotted as dotted lines (ma- 
genta) and large amorphous grains as dashed lines (magenta). 



approximation. For EP Lyr the GRF grains prove the best match, 
but the difference in x 2 with the small DHS grain approxima- 
tion is only minimal. The best fit to EPLyr is obtained using 
both small (0. 1 /mi) and larger (2.0 /mi) silicate grains. The small 
difference in calculated x 2 values for EP Lyr is due to the low 
signal-to-noise ratio, making it hard to distinguish between dif- 
ferent synthetic emission profiles. For HD 52961 small grains in 
Mg-rich silicates give the best^- 2 . Plots of our best fitting models 
can be found in Fig. [2] Table[3]gives the resulting parameters. 

The large x 2 value of HD 52961 quantifies that this star has 
a very peculiar, unique chemistry, and we did not succeed in 
explaining all of the observed features. The strong forsterite 
11. 3 /mi feature in the GRF dust approximation is clearly too 
broad. DHS grains fit the feature better, but other feature profiles 
are fitted less well with this approximation. There also appears 
to be a short wavelength shoulder on the amorphous 9.8 /im fea- 
ture, which is not explained in the modelling. The strong 16.5 /mi 
feature is not reproduced in central wavelength by any of the 
different m odels. We already o bserved this trend in our full sam- 
ple fitting dGielen et al. 2008), where the feature seemed to be 
shifted bluewards in comparison with the mean spectrum of the 
full sample. The two narrow features around 19//m could be an 
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Fig. 3. The Spitzer-IRS spectrum of EPLyr compared to the 
ISO-SWS spectrum of IRAS 09425-6040. The spectrum of 
IRAS 09425-6040 is normalised and offset for comparison. 

Table 2. x 1 values for different models used in our full spectral 
fitting. 





EPLyr 

x 2 


HD 52961 

x 2 


model description 






model 1 


21.7 


129.4 


Mie-0.1 -2.0jum- 


A" = 


= 0.5 


model2 


6.2 


67.5 


DHS-0.1 - 1.5 /an 


- X 


= 0.5 


modeB 


6.2 


63.8 


DHS-0.1 - 1.5 yum 


- X 


= 1.0 


model4 


8.4 


101.4 


DHS - 1.5- 3.0 ^m 


- X 


= 0.5 


model5 


8.6 


140.8 


DHS - 1.5- 3.0 ^m 


- X 


= 1.0 


model6 


5.9 


64.2 


GRF-0.1 -2.0yum 


- X 


= 0.5 


model7 


6.3 


50.0 


GRF-0.1 -2.0yum 


- X 


= 1.0 


model8 


5.8 


96.5 


GRF-2.0-4.0yum 


- X 


= 0.5 


model9 


5.4 


72.2 


GRF-2.0-4.0/um 


- X 


= 1.0 



Note: For each model we give the used dust approximation, grain size 
and Mg-Fe content in the amorphous grains, x = 1.0 denotes pure 
Mg-rich amorphous dust, x = 0.5 the more standard Mg-Fe amorphous 
silicates. 



artifact of the data reduction, since in this region there can be a 
bad overlap between the SH and LH Spitzer-IRS high-resolution 
bands. 

For EPLyr we fit the spectrum longwards of 18 fim, where 
the silicate features are seen. This gives dust temperatures be- 
tween 100 and 230 K. This model, however, does not fit the spec- 
trum before 18 /mi, since the continuum does not follow the ob- 
served strong downward slope before 20 /mi. If we try to fit the 
full Spitzer wavelength range we find we can get a better fit to 
the underlying continuum but then the features at 27 and 33 /mi 
are much stronger in the observed spectrum than in our best 
model. Unlike in other sources, a two temperature approach fails 
to model both the observed continuum and the coolest features 
for the full Spitzer wavelength spectrum of EP Lyr. Irrespective 
of the derived continuum temperature, all the tested models 
give estimates of the dust temperatures between 100 - 300 K, 
which agrees with the temperatures derived in the SED mod- 
elling (Sect. [3J. Clearly, the crystalline dust particles must be 
quite cold. 




o.o h ... i ... i ... i ... i ... i ... i ... -I 

4 6 8 10 12 14 16 18 

Wavelength [jj.m] 



Fig. 4. Comparison between the continuum-subtracted spectrum 
of EPLyr and HR4049, another mixed chemistry source. Both 
stars show clear emission due to CO2 gas around 15//m and 
strong PAH features before 13 /mi. The CO2 features at 14.9 - 
16.2 jum, 15.3 /mi and 15.1 /mi are respectively due to 12 C 16 C»2, 
13 C 16 2 and 16 12 C 18 0. The PAH emission is clearly very dif- 
ferent in HR4049, where class B PAHs are found, whereas 
in EPLyr the PAH features can be attributed to class C (see 
Sect. 14741 



4.3. C0 2 emission 

4.3.1. Introduction 

CO2 emission has been found in approximately 30% of all 
O-rich AGB sta rs dJusttanont et all I1998I: iRvde et all 1 19991 
ISloan et al.l2003l) . but CO2 detections in post-AGB stars are rare. 
To our knowledge CO2 gas has been foun d in only two post- 
AGB stars, the Red Recta ngle and HR4049 d Waters et alJI 199 81: 
ICami & Yamam ura 2001), which are also binaries surrounded 
by a stable circumstellar disc. HR 4049 is the only example of a 
post- AGB star showing CQ 2 in emission in the 13 - 16 /mi re- 
gion. Cami & Yamamura] d2001l) argued that the isotopic distri- 
bution of oxygen in HR 4049 is abnormal, based on the isotope 
ratio anal ysis of the CO2 emi ssion features. This was not con- 
firmed by iHinkle et al.l {2007), who use high-resolution spectra 
of the fundamental and first overtone CO vibro-rotational transi- 
tion in the near-IR. 

Both EPLyr and HD 52961 show clear gas phase emission 
lines of 12 CC>2 and 13 CC<2. These emission li nes were also seen i n 
only one other source in our Spitzer sample (Gi elen et al.l l2008). 
namely in IRAS 1 1 74-5704. 

The CO2 emission of EP Lyr seems to be lying on top of a 
"plateau" that extends from 13 to 17 /mi. A similar plateau i n this 
region is observed in PAH-rich sources dPeeters et al. 2004), but 
this plateau is much broader and ranges from 15 to 20 /mi and is 
often characterised by strong emission features at 16.4/mi (and 
less prominent at 15.8, 17.4 and 19/mi), and thus quite different 
to the one seen in EP Lyr. 

4.3.2. Analysis 

To retrieve the very rich spectral information of the CO2 emis- 
sion bands, we calculate s pectra of CO2, using HITRAN line 
lists (Rothman et al. 20 051) and a circ ular slab model for the 
radiative-transfer dMatsuura e t al. 2002). The model has four pa- 
rameters: the excitation temperature (T ex ), the total CO2 col- 
umn density (AO, the radius of CO2 gas (r) and the isotope ra- 
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Table 3. Listed are the best fit parameters deduced from our full spectral fitting. 



Name 


x l 




Tdustl 


Fraction 


Tcontl 


Tcontl 


Fraction 






(K) 


(K) 




(K) 


(K) 


T cont 1 ~ Tcontl 


EPLyr 
HD 52961 


5.4 
50.0 


200^° 


228 48 " 
774186 

' ^96 


O9O.I _Ql0.6 
O9 8.0_ o [8.1 


205 702 
111 358 


641 331 

u -301 

996 4 


96° -02 - 04 UU6 

99O.OO _ O1 8.05 



Name 


Olivine 


Pyroxene 


Forsterite 


Enstatite 


Continuum 




Small - Large 


Small - Large 


Small - Large 


Small - Large 




EPLyr 


6 41 - 5 2i 






5f - 28^ 




HD 52961 


13 - l 41 


55 12 -2 43 


6^-33- 


l 8 - 3 26 ~ 


69 4 



Note: Top part: The^ 2 , dust and continuum temperatures and their relative fractions. Bottom part: The abundances of small (0.1 yum) and large 
(2.0yum) grains of the various dust species are given as fractions of the total mass, excluding the dust responsible for the continuum emission. The 
last column gives the continuum flux contribution, listed as a percentage of the total integrated flux over the full wavelength range . The errors 
were o btained using a Monte-Carlo simulation based on 100 equivalent spectra. Details of the modelling method are explained in Giel en et"aT] 
J20081) . 



tio. The radius of the CO2 layer is given relative to the radius 
of the background continuum source at 13/um. The dependence 
of the CO2 model spectra on these parameters are described by 
Cami & Yamamura (2001). We estimate a pseudo-continuum by 
using a spline fit and a linear fit for HD 52961 and EPLyr, re- 
spectively. For HD 52961, we interpolate the spectrum at the 
spectral range where CO2 bands have little influence on the ob- 
served spectrum. The continuum was also chosen so that the 
forsterite feature at 16^/m would be removed. A spline fit was 
tested for EP Lyr but failed because of the richness of CO2 fea- 
tures in the mid-infrared range, so we simply use a linear in- 
terpolation between 13.4 and 17.8/mi. Estimated continua are 
displayed as dotted lines in the top panels of Figs. and [6] The 
resulted parameters for the model calculations are summarised 
in Table |4] and the resulted spectra show the identifications of 
the different CO2 bands (bottom panels of Figs. [5] and |6j . 

Many small features in EPLyr in the 13.5 - 17//m region 
are due to C0 2 : features at 13.5, 13.9, 14.7, 14.9, 16.2 /zm are 
attributed to the CO2 main isotope I2 C 16 C>2. The main isotopic 
12 C 16 C>2 bands are probably optically thick, surpressing the line 
intensities. 13 C 16 2 and 16 12 C 18 bands are found at 15.3/im 
and 15.1 yum, respectively. The prominent I6 12 C I8 feature is 
surprising. This fea ture was also found in the other binary post- 
AGB star HR 4049 (|CamL&^mamr3lQQl. 

The model uses a high fraction of isotopes, but actual abun- 
dance ratios remain largely uncertain, mainly because of the un- 
certainty of the interpolated continuum spectrum and the optical 
thickness of the main isotope. Nevertheless, these two features 
are particularly prominent in the spectrum of EP Lyr, more than 
in HD 52961, suggesting different isotope ratios for EPLyr. We 
see that the observed "plateau" in EP Lyr can be explained by the 
richness of the 12 C 16 C>2 features, but the 12 C7 I3 C ratio is con- 
firmed to be low. The strength of the isotopes, including the very 
rare ls O (in the Sun 16 0/ 18 is ~ 500), is an exclusive feature of 
post-AGB stars. 
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Fig. 5. EPLyr: top panel shows the observed spectrum (solid 
line) and pseudo-continuum spectra (dotted line). The bottom 
panel shows the continuum divided spectrum, CO2 model spec- 
tra (combining all of the isotopes) and individual CO2 isotope 
spectra (from top to bottom 12 C 16 2 , °C 16 2 , 16 12 C 18 0). 

Table 4. Resulting parameters for the model calculations. The 
isotope ratio is as follows: 12 C I6 2 : 13 C 16 2 : 16 12 C I8 0. 



4.4. PAH features 

Polycyclic aromatic hydrocarbons are found in a large variety 
of objects, including the diffuse ISM, HII regions, young stellar 
objects, post-AGB stars and planetary nebula e. They have s trong 
emission features in the 3-13^m region (e.g. Tielens 2008). The 
feature at 3.3 yum is arises from the C-H stretching mode of neu- 



rit N r isotope ratio 

(K) (cm' 2 ) (R,) 

EPLyr 900 8 x 10 1!i 4?7 0.7 : 0.2 : 0.1 

HD 52961 800 5 x 10 18 4.7 0.93:0.05:0.02 
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Fig. 6. Same as Fig. |5]but for HD 5296 1 . 



tral PAHs. The C-C modes produce features with typical central 
wavelengths at 6.2 and 7.7 jjm. The 8.6 /mi feature is due to C-H 
in-plane bending modes and features longward of 10 /im can be 
attri buted to C-H out-of -plane bending modes. 

iPeeters et al.l (120021) defined three groups of PAH spectra 
based on their emission profiles and peak positions. The "class 
A" sources have features at 6.22, 7.6 and 8.6 fim. "Class B" 
sources show the same features but shifted to the red, peaking 
at 6.27, 7.8 and > 8.6 fim. They also identified two "class C" 
sources, the Egg Nebula (AFGL2688) and IRAS 13416-6243, 
both post-AGB objects. These rare "class C" sources show emis- 
sion features at 6.3 jum, no emission near 7.6/im, and a broad 
feature centred around 8.2 fim, extending beyond 9 fim. 

With the release of the IRS aboard the Spitzer Space 
Telescope, a limited number of additional class-C sources were 
discovered (Fig. |7]>. MSX S MC 029, a class-C post-AGB star 
in the SMC, was detected bv lKraemer et al.l d2006l) . ISloan et alj 
(2007i) report on the detection of class-C PAH features in 
HD 100764, a carbon-rich red g iant with evidence of a cir- 
cumstellar disc. lJura et all (J2006) also report on the detection 
of class-C PAH features in a circumstellar disc around the 
oxygen-rich K-giant HD 2335 17 . Two young objects, the T Tauri 
star SUAur (Furla net al.l 120061) and the Herbig Ae/Be source 
HD 135344 dSloan et alj|2005l) . also show PAH spectra of class 
C, although in HD 135344 the PAH features seem to be some- 
what more in between B and C. This source is also slightly hotter 
than other class C sources. A comparison of the PAH features in 
all these sources is given by Slo an et all d2007l) . They find that 
all the known class-C spectra are excited by relatively cool stars 
of spectral type F or later and argue that the hydrocarbons in 
these sources have not been exposed to much ultraviolet radia- 
tion. The class-C PAHs are then relatively protected and unpro- 
cessed, while class A and B PAHs have been exposed to more 
energetic photons and are hence more processed. 




6 8 10 12 14 

Wavelength [fir 

Fig. 7. Continuum- subtracted s pectra (based on a spline fit) of 
class-C sources as described in Slo an et al. (2007). The vertical 
dashed lines are at 8.1, 11.3 and 12.6/mi. For comparison we 
also plot HD 44179, which has class B PAH emission. 
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Fig. 8. Comparison between the Spitzer-IRS spectrum of 
HD 233517 and EPLyr in the 6-14^111 region. 



EPLyr shows PAH emission bands at 8.1, 11.3 and 12.6//m. 
Other features can be seen at 13.25jt/m. Whether these can be 
attributed to PAHs remains uncertain. The PAH spectrum of 
EP Lyr can be classified as class C. EP Lyr is a high amplitude- 
variable with an effective temperature around 7000 K, which is 
on th e hot end of the other class-C emitters (see ISloan et ail 
12007b . 

Sofar, it is not clear whether the PAH carriers reside 
in the circumbinary disc, or in bipolar lobes created by 
a more recent mas s- loss event, as obse r ved in HR4049 
19991: iDominik et all 120031: lAntoniucci etaTJ 



120051: iHinkle et all 120071: iMetiut et alJ 120091) and in the Red 
Rectangle dMen'shchikoy et all 120021: ICohen et all 12004). For 
HR4049. Hohnson et all d 19991) found that the optical polarisa- 
tion seems to vary with orbital phase, suggesting that the polari- 
sation in the optical is due to scattering in the circumbinary disc. 
In the UV, the polarisation is caused by scattering in the bipolar 
lobes, which should contain a population of small grains, includ- 
ing the PAH carriers. HR4049 and HD 44179, as well as EPLyr 
and HD 52961, are strongly depleted, and molecules or the for- 
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mation of dust in these very depleted environments is likely very 
different from solar condenstation. As the CO molecule is abun- 
dant in the circumstellar environment, accretion of circumstel- 
lar gas will likely result in a C~0 photosphere. HR4049 and 
HD 44179 are stars with PAH emission belonging to the more 
standard class B. 

Unlike what is detected around HD44179 (the ERE neb- 
ula) and HR 4049 (nano-diamonds), for EP Lyr the PAHs are the 
only carbon-rich component observed in the circumstellar spec- 
tra. There is no evidence that the photosphere of EP Lyr is in, or 
went through, a carbon-rich phase. The photosphere is depleted 
so that nucleosynthetic yields are very hard to recognise, bu t 
the photosphere is clearly oxygen-rich dGonzalez et al.|[T997l) . 
A scenario involving hot-bottom burning to return to an oxygen- 
rich condition after a carbon-rich phase in the stellar evolution 
along the AGB is very unlikely: it would imply the object is of 
a mor e massive origin dLattanzio et alJ[T996t iMcSavenev et alj 
12007b . which is in contradiction with its Galactic coordinates 
of b = 6.9°. Assuming 10000L G for a putative massive pro- 
genitor, the distance above the Galactic plane would be about 
660 pc, which is very high for an object of 5 - 6 M Q . Moreover, 
the initial metallicity is likely subsolar as indicated by the solfur 
and zinc abundances. Although, with the depletion, it is unclear 
whether these abunances of the volatiles do indeed represent the 
initial conditions. Also the pulsation period of EP Lyr is similar 
to other RV Tauri objects which are of low initial mass. We con- 
clude that the photosphere of EP Lyr is now O-rich, and we argue 
it is very unlikely that it was ever in its history in a carbon-rich 
phase. The PAH synthesis likely occurred in O-rich conditions. 

As Figure [8] shows, the spectrum of EP Lyr has a s triking re- 
semb la nce t o that o f HD 233517, shortward of 13 fim Oura et alj 
2006). Euri d2003l) hypothesises that HD 2335 17 was a short- 
period binary on the main sequence. A circumstellar disc was 
then formed when the companion star was engulfed by the more 
massive star when it entered giant evolution, followed by a phase 
of strong mass ejection in the equatorial plane. Since HD 2335 17 
is an oxygen-rich star, it is re markable that the disc shows fea- 
tures of carbon-rich chemistry. lJura et all ((2006) propose a sce- 
nario in which the PAHs could be formed inside the disc due 
to Fischer- Tropsch (FT) catalysis on the surface of solid iron 
grains. These FT reactions ca n convert CO and H2 into water and 
hydrocarbons dWilla cv 2004), these hydrocarbons could then be 
converted in into PAHs. 

So far it remains unclear whether the shape of the observed 
emission features can detect if the PAH carriers reside in the disc 
or an outflow. 

5. Spectral energy distribution 

5.1. 2D disc modelling 

For both stars extensive photometric data are available. This, to- 
gether with the Spitzer infrared spectral information, allows us 
to constrain some of the physical characteristics of the circumbi- 
nary disc. 

We fit the SED using a Monte Carlo cod e, assuming 2D- 
radiat i ve-transfer in a passive disc model dDullemon d et al.l 
I2OOU iDullemond & Dominikl |2004|) . This code computes the 
temperature structure and density of the disc. The vertical scale 
height of the disc is computed by an iteration process, demand- 
ing vertical hydrostatic equilibrium. The distribution of dust 
grain properties is fully homogeneous and, although this model 
can reproduce the SED, dust settling timescales indicate that set- 
tling of large grains to the midplane occurs. Using the dust set- 



tling timescale 

n E 1 z 

t se t = x FT ln — 

2 p d a Ll k zo 

with Zo the surface density, p,/ the particle density, a the 

grain size and Q.k = the Keplerian rotation rate 

dMivake & Nakagawal H995), we find that grains with sizes be- 
tween 500 fim and 0. 1 cm can descend 50 AU on timescales sim- 
ilar to the estimated lifetime of the disc. An inhomogeneous disc 
model wit h a vertical gradien t in grain-size distribution is thus 
necessary (Gi elen et al.l l2007). These large and cooler grains in 
the disc midplane mainly contribute to the far-IR part of the SED 
and constitute the main fraction of the total dust mass. The disc 
structure and near- and mid-IR flux are almost fully determined 
by small grains. So we use a homogeneous 2D disc model to fit 
the near- and mid-IR part of the SED and add a single black- 
body temperature to represent the cooler midplane made up of 
large grains. The lack of observational constraints on the temper- 
ature structure of this component does not allow us to constrain 
this extra parameter. 

Stellar input parameters of the model are the luminosity, 
the total mass (we assume the total gravitational potential to be 
spherically symmetric with a total mass of M = 1 M ), and T e ff . 
The luminosity (and thus the distance) for these sources is not 
well constrained so we use values between L = 1000 - 5000 L , 
typical values for post-AGB sources. Input disc parameters are 
R m and R out , the different dust components, the total disc mass 
and the power law for the surface-density distribution. Since we 
are not dealing with outflow sources a power law a > -2 is used. 
The gas-to-dust ratio is kept fixed at 100. 

The modelling is still degenerate, especially in parameters 
like the outer radius and the total disc mass which can be eas- 
ily interchanged, without strongly influencing the SED. We use 
a dust mixture of amorphous and crystalline silicates in grain 
sizes ranging of 0.1 - 20 fim, with a power law distribution of 
-3.0, for the homogeneous disc. For HD 52961 the 850 fim sub- 
millimetre data points to the presence of extremely large grains 
in the disc, but these large grains are assumed not to influence the 
near- and mid-IR part of the SED, and will only be important for 
the blackbody component. We use a value of 300 AU for R out . 
For the inner rim we use the radius at which the temperature of 
the inner rim equals about 1500 K. This is a typical value for the 
dust sublimation temperature of silicates, although values as low 
as 1 200 K are sometimes also used. The total SED energetics are 
then calculated, given a specific inclination angle of the system. 

We do not aim to reproduce the observed narrow features in 
the Spitzer spectrum, since these features originate from an op- 
tically thin upper layer of small grains at the disc surface. They 
can only be fitted well using an inhomogeneous disc model with 
grain settling. Instead, we want to model the observed general 
energetics of the SED spectrum, thus the observed continuum 
and amorphous dust features. 

5.2. Comparison with interferometric data 

The circumstellar environment of HD 52961 has been resolved 
using the VLTI/MIDI instrument, with angular sizes in the N- 
band between 35 mas a nd 55 mas in a uniform disc approxima- 
tion dDeroo et alj|2006h . EP Lyr is too faint for current interfero- 
metric capabilities. 

To compare our disc model of HD 5296 1 with the MIDI data 
we made model images from which we extracted visibilities, us- 
ing the same projected baseline lengths (40m/46m) and angles 
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Fig. 9. Some SED disc models showing the influence of the in- 
clusion of metallic iron and different grain sizes in the fitting. 
All models have the same physical parameters, the only differ- 
ence being the amount of iron in the disc and the grain sizes. The 
solid line depicts a pure silicate disc with grain sizes between 
0.1-20yum. The dashed line represents a model with 5% metallic 
iron and 95% silicate with grain sizes between 0.1-20^m. The 
dotted line represents a silicate disc model with grain sizes be- 
tween 0.1-50/mi and the dot-dashed line one with grain sizes 
between 0. 1 - 1 00 /mi. 



(45°/46°) as the observations. The only free parameters here are 
the inclination of the disc and the orientation angle of the system 
on the sky. A range of inclinations which still fit the observed 
SED was tested, in steps of 15°. The orientation angle is var- 
ied 1° at a time. The result of this comparison can be found in 
Sect. 13311 

5.3. Results 

When modelling the near- and mid-IR part of the SEDs, the 
feature-to-continuum ratio of the silicate features is too strong in 
comparison with the observed spectra for both stars. Moreover, 
th e near-IR flux is of ten underestimated. This was also observed 
in lGielen et al.l d2007l) . where we fitted the SED similarly to two 
post-AGB stars, RU Cen and AC Her. Including an extra contin- 
uum opacity source is needed to reduce the strength of the fea- 
tures (see Fig. [9]» a nd to increase the n ear-IR contribution. From 
our previous work dGielen et al.l200 8) and Section l4~2l we found 
evidence that (a fraction of) the silicates might be iron poor. So 
we use metallic iron as a potential opacity source: its near-IR 
opacity is large, but the absorption coefficient is unfortunately 
featureless so direct detection is difficult. Inclusion of free metal- 
lic iron has a strong impact on the modelling because the near-IR 
excess increases significantly with a given inner radius. Another 
possible opacity source is the inclusion of hot large grains in the 
homogeneous disc model. 

5.3.1. HD 52961 

We can use the interferometric data for HD 52961 to constrain 
the distance to the system (Fig.lTTb. When we compare the mod- 
elled visibilities with the observations for HD 52961 we find that 
the visibilities of our model are too high, when using standard 
models of L — 5000 L which fit the SED and impose the in- 
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Table 5. Results of our SED disc modelling. 
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Note: Given are the inner and outer radius (Ri n -R ut), the total disc mass 
m for the homogeneous disc model, the surface-density distribution 
power law a, the percentage of iron in the homogeneous disc model, 
the blackbody temperature and the inclination of the system. 



ner radius to be at sublimation radius. This means we need to 
increase the angular size of the N-band emission region, either 
by increasing the physical size or by decreasing the distance. 
Increasing the disc size to outer radii > 500 AU does not in- 
fluence the N-band emission so we need to increase the inner 
radius, flatten the density distribution power law and/or decrease 
the luminosity in the disc model. Changing the surface-density 
distribution power law to values > —1.5 proved incompatible 
with the observed SED. We therefore use an average luminos- 
ity of 3000 L Q and move the inner radius to larger distances. A 
good fit to the SED was obtained using an inner radius of 10 AU. 
This assumed luminosity gives a distance to the system of about 
1700pc. At 10 AU the temperature of the inner rim is around 
HOOK, which is slightly below the canonical dust sublimation 
temperature for silicates. 

For HD 52961 we tested models with and without the inclu- 
sion of metallic iron. The resulting fits can be seen in Figure [TOl 
parameters can be found in Table|5] Since this star shows a rather 
strong 10 /.im feature, pointing to relatively large amounts of hot 
dust in the disc, we find we need a rather steep surface-density 
distribution, a < -1. Since we only have one flux point at long 
wavelengths, we add a simple 160 K blackbody model to fit the 
observed submillimetre emission. 

When no iron is included (model A) we see that the flux 
around from 2 to 8 fim is strongly underestimated. Including 
about 10% metallic iron (model B) increased this flux signifi- 
cantly. The inclusion of metallic iron has only a minor influence 
on the N-band interferometric measurements. It decreases the 
modelled visibilities by about 10%, which is still consistent with 
the observed visibilities. 

The modelled visibilities (Fig. ITTb lie within the observed 
visibility range but a remarkable detection is that the variation in 
visibilities between the two observations is quite large, despite 
the very limited difference in lengths (41.3- 46.3 m), as well as 
in projected angle (45.6-46.3°). This is illustrated when we plot 
the visibility versus the spatial frequency for a given wavelength, 
as seen in Figure [12] In this figure we illustrate that when using 
a uniform disc for the intensity distribution, the steep visibility 
drop can be accounted for. The physical disc model is, however, 
much smoother and does not contain the very sharp edge charac- 
teristic of the uniform disc. The model also do not rep roduce the 
observ ed 'bump' in visibility between 9 and 12 fjm. iDeroo et al.l 
(2006) explain this observed increase in visibility as being due 
to a non-homogeneous distribution of the silicates, which con- 
tribute most to the inner regions of the disc. The current disc 
model does not include the physics to be able to reproduce this 
radial distribution of dust species. 



10 



C. Gielen et al.: Analysis of the infrared spectra of the peculiar post-AGB stars EPLyr and HD 52961. 



<u 10 




000.0 



<u 10 




10.0 100.0 
A [yum] 



1000.0 



Fig. 10. SED disc modelling of HD 52961 (top: model A without 
metallic iron, bottom: model B with metallic iron). The dashed 
line represents the homogeneous disc model consisting of grains 
between 0. 1 fim and 20 /urn. The solid line gives the disc model 
with an added blackbody to represent the cool midplane. Crosses 
represent photometric data and in the infrared we overplot the 
observed Spitzer-IRS spectrum. The dotted line represents the 
adopted Kurucz stellar model. 



The submillimetre 850 //m flux for HD 5296 1 and the derived 
blackbody temperature of 160 K can be used to estimate the dust 
mass of large grains in the disc. In t he optically thin approach 
the disc mass can be estimated using (iHildebrandl 19831) 



M rl = 



^850 D 

k$50 Bg 5 o(T) 



Assuming a cross section of large spherical grains, the mass 
absorption coefficient of 850yum grains in blackbody approx- 
imation is about 2.4 cm 2 g -1 . The mass absorption coefficient 
is given by k — 4™ 3 , with a the grain size, and p typically 

3.3 gem -3 for astronomical silicates. This results in dust-mass 
estimates in large grains of 3 .2 x 1 -6 M for HD 5296 1 . 



5.3.2. EPLyr 

The SED-fitting gives an estimate of the distance to the system 
of d = 3.4 kpc, assuming a luminosity L = 3000 L . This estima- 
tion is not only dependent on the assumed luminosity of the star, 
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Fig. 11. The observed MIDI measurements of HD 52961 (gray 
data with error bars) and the visibilities deduced from disc model 
A (black solid line), for the two different baselines. The model 
has an inclination of 45° and a position angle on the sky of 227° 
East of North. 
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Fig. 12. Comparison between modelled and observed visibilities 
at 12.6 /zm, at the two different baseline lengths. The dashed line 
represents a uniform disc model with an angular size of 50 mas. 
The solid line gives the visibilities as calculated from the 2D disc 
modelling (model B). 



but also on the adopted inclination of the system. Other derived 
disc parameters can be found in Table [5] 

For EPLyr it proved very problematic to get a good fit 
to the observed strong 20 fim feature, without introducing a 
strong 10/im silicate feature. The mixed chemistry adds to 
the complexity of the object as the PAH emission and under- 
lying continuum in the near-IR may very well come from a 
distinct structural component, in for example the polar direc- 
tion. This is already seen in HD 44179, w ere the observed PAH 
emission comes from a bipolar outflo w (Breg man et al.l 1 1993; 
iMen'shchikov et alj|2002t ICohen et alj|2004l) . The lack of data 
shortwards of 7 /jm makes it hard to get a good continuum esti- 
mate of the near- and mid-IR energetics. Cool dust clearly dom- 
inates the SED, but without additional information from inter- 
ferometry, we cannot constrain parameters like the inner radius 
or the surface-density distribution. We opted to keep a rather flat 
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Fig. 13. The SED disc model of EP Lyr. The model parameters 
are given in Table [5] 



density distribution of a = -1.0 and do not force the model to 
fit the spectrum shortward of 15 fim. 

To get a good fit to the SED a very large inner radius of 
about 40 AU is necessary. At this distance the inner rim reaches 
temperatures of ~ 300 K. This temperature agrees with the tem- 
peratures found in the spectral modelling (Sect, \4.2i . If we in- 
clude metallic iron in the model we find we need inner radii 
even larger than 200 AU, which seems physically implausible. 
Unfortunately we do not posses submillimetre data for this star 
so we cannot determine the blackbody temperature associated 
with the midplane. The small Lir/L* as well as the lack of a 
near-IR dust excess shows that the inner rim is quite far from the 
sublimation radius. 



6. Conclusions 

HD 5296 1 and EP Lyr both have rich infrared spectra, and the as- 
sembled multi-wavelength data show that these evolved objects 
are surrounded by a stable circumstellar disc. While the binary 
nature of HD 52961 is well established, the binarity of EPLyr 
is suspected but not yet firmly proven. The discs must be cir- 
cumbinary as the sublimation radius of the dust is larger than the 
determined (HD 52961) and suspected (EP Lyr) binary orbit. 
Recent studies have shown that many of these binary 



y 

post-AGB systems ar e already detected (IVan W inckel 2003; 
De Ruvter etail l2006t iDeroo et all l2006l l2007at iGielen et all 



2008), but EPLyr and HD 52961 both show quite distinct char- 
acteristics in dust and gas chemistry as well as in physical prop- 
erties of their discs. 

EPLyr and HD 52961 are the only stars from the larger 
Spitzer sample that have clear CO2 gas emission lines in the mid- 
IR. Our modelling shows that the emission in both stars can be 
well fitted and is dominated by 12 C I6 C>2 features, but clear de- 
tections of other isotopes are present as well. Similar excitation 
temperatures and column densities are found in both objects, but 
with different ratios for 12 C 16 2 and 13 C 16 2 . Why these two 
stars are the only ones from the larger sample showing strong 
CO2 features, and if there is any relation with the low observed 
infrared flux remains unclear. Similar feature strengths observed 
in the other stars would have been easily detected. One possibil- 
ity is that the low dust emission in the two sample stars, reflected 
in the low LjB./L star , makes it easier for the CO2 gas to become 
visible. This effect is also seen in AGB stars, where CO2 gas 



emission is strongest in sources with the lowest mass-loss rates 
(ISloan et alJll995 10^1200% ISloan et aTJl2003T) . 

One of the most remarkable features is the clear detection 
of 18 Q isotopes of CQ 2 in both objects. Together with HR 4049 
(ICami & Yamamurall200ll) . the strong 16 12 C 18 band is a sys- 
tematic feature of the gas emission in the discs of post-AGB 
binaries when CO2 emission is detected. The high 18 abun- 
dance of HR4049 derived in an opti cally thin approxim ation 
of a putative nucleosynthetic origin (Lugaro_etaL|[2005) was 
not confirmed by t he analysis of CO first overtone absorption 
dHinkle et al.ll2.007h in the same object. It is likely that the CO2 
gas is strongly optically thick, also in EPLyr and HD 52961, so 
that very rare isotopes are detectable. 

The high-resolution Spitzer spectra also reveal unique solid 
state features. As in the bulk of the disc sources (Giele rTet al.l 
120081) . crystalline silicate features prevail in both stars, but unlike 
what we found for the larger sample, they proved very hard to 
model. In HD 52961 we observe some unique strong crystalline 
features at 11.3 and 16 /mi, which could not be reproduced in 
the modelling, irrespective of the grain size used in the mod- 
els, shape or assumed grain model. In the 2D disc modelling we 
could not fit the steep rise around 10 /im, without the inclusion of 
metallic iron. Combining our physical model constrained by the 
SED, together with our interferometric data, we concluded that 
the inner dust rim is slightly beyond the dust sublimation radius. 
This is in contrast to similar binary objects like IRAS 08544- 
4431 where the interferometric data shows that th e dusty disc 
has to start very near to the dust sublimation radius (IDeroo et al.l 
l2007al) . Assuming a luminosity of 3000 L , we find that the in- 
ner disc radius of HD 52961 is rather large, around 10 AU. The 
strong 850/mi flux shows that this object has a component of 
very large grains. This contribution was added to the SED fitting 
by an additional colder Planck curve. 

EP Lyr has only a very small infrared excess, but the Spitzer 
spectrum is very rich in spectral details. The most remarkable 
characteristic is the clear PAH emission, in combination with 
the strong crystalline features at longer wavelengths. There is no 
evidence that the central star evolved into a carbon star when on 
the AGB, yet unprocessed class-C PAH features are clearly de- 
tected. Whether these PAH species are formed in the circumbi- 
nary disc or in a recent, likely polar outflow of the depleted cen- 
tral star, remains unclear. An extra component of cold dust is 
necessary in this object as well, to fit the entire Spitzer spectrum. 
Unfortunately, EP Lyr is too faint for the current interferometric 
possibilities. 

The mixed chemistry, the strongly processed cold crystalline 
silicates and low 12 C/ 13 C ratio are in common with the subgroup 
of silicate J-type carbon stars, which can also display strong 
crystalline material. This corroborates the conclusion that in the 
latter, the disc is circumbinary. The abundance studies of J-type 
carbon stars are not complete enough to probe whether photo- 
spheric depletion affected these objects as well. 

Both objects are extreme examples of post-AGB binary stars, 
with characteristics dominated by the presence of a stable cir- 
cumbinary disc. This disc environment is, to first order, well 
modelled by assuming a passive, irradiated stable disc. In this 
paper we corroborate that this geometry is ideal to induce strong 
grain processing and a rich, even mixed chemistry. We conclude 
also that a homogeneous disc model is too primitive to model the 
spectral details as evidence of grain settling is strong. The route 
to PAH formation (and excitation) in the O-rich EP Lyr remains 
to be studied in detail as PAH emission is only observed in a very 
limited number of such sources. 
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This detailed study of HD 52961 and EPLyr shows that 
many questions still remain in our current understanding of the 
evolution of a significant number of post-AGB binary stars, and 
the impact of the circumbinary discs on the entire system. 
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